We describe a linear Bayesian inversion method to estimate the relevant petrophysical properties of the media forming a reflecting interface from the observations of amplitude variation with incidence angle. Three main steps characterize the proposed approach: – information from borehole logs are statistically analysed to estimate the empirical models that describe the functional relationship between petrophysical (e.g. porosity, saturation, pressure or depth) and seismic variable(P and S velocities and density); – the pure-mode (PP) reflection coefficient is parameterized in terms of the relevant petrophysical variables and is linearized in order to implement the linear inversion; – the sought petrophysical parameters are estimated from the seismic reflected amplitudes by applying the linearized inversion where a priori information, data and model errors and solutions are described by probability density functions. We test the method on synthetic and real data relative to reflections from a shale/gas-sand interface where the amplitude versus angle response, besides the lithological contrast, is mainly controlled by the saturation and porosity of the sand layer. The outcomes of the linearized inversion are almost identical to those obtained by a previously developed non-linear inversion method demonstrating the applicability of the linear inversion. It turns out that the gas-sand saturation in the range 0%–95% is a poorly resolved parameter while the porosity is the best resolved parameter. The issues of robustness and resolution of the inversion are discussed either through singular value decomposition analysis or the observation of the a posteriori probability density functions. The linear inversion algorithm, compared with the previously developed non-linear method, reduces significantly the computation time allowing for more extensive applications. 相似文献
Side scan sonar records, sediment textural characteristics, and in-situ field observations were used to study gravelly and pebbly dunes and sorted bedforms on the inner shelf of Marettimo Island, along the northwestern Sicilian shelf. The dunes are composed of coarse sands, gravels and pebbles (D50: 2–16 mm), displaying a symmetrical shape with a wavelength in the range of 1–2.5 m and a height of 0.15–0.30 m. The bedforms are distributed in a patchy pattern in a depth range of 10–50 m, and are described for the first time on a Mediterranean inner shelf. Sorted bedforms are linear morphological features developed almost perpendicular to the coast in the eastern sector of the island between 15 and 50 m water depth. Bottom shear stresses required for sediment entrainment and the generation of the shallower dunes can be reached during strong storms (Hs=5–6 m; Tp=9–11 s), which are not common in the Mediterranean Sea. However, wave storm events recorded in the study area during the last 17 years are not able to generate the coarsest and deeper dunes, suggesting that the stirring mechanism for dune formation is associated with severe storms that have a recurrence interval of more than 17 years. The long-term stability of the coarse bedforms is supported by the permanence of sorted bedforms without significant morphological changes for long periods (>13 years). Therefore, it is shown that processes forming coarse bedforms can occur in tideless and moderate-energy settings like those of the Mediterranean continental shelves, although the morphological features are probably less dynamic and remain unaltered for longer periods than on higher-energy shelves. 相似文献
The study proposes a model by which a thick succession of volcanic tuffs can be zeolitized by alteration of pyroclastic material
in the presence of sufficient eruptive water and at temperatures close to water vapour condensation. In the case of phreatomagmatic
products, the model simplifies interpretation of problematic deposits that exhibit pronounced vertical and lateral variation
in lithification grade. A major feature of the model is that thick zeolitized tuffs can be formed during emplacement of pyroclastic
products, in marked contrast to later alteration in an open hydrologic system. Geological, volcanological and mineralogical
data for the Neapolitan Yellow Tuff, a widespread trachytic pyroclastic deposit outcropping around Campi Flegrei (Southern
Italy), have been used to infer the physico-chemical conditions that determined mineral genesis. This tuff shows a reduction
in lithification grade towards the base, top and with distance from the vent and very variable zeolitization within the lithified
portion. We suggest that during initial emplacement the erupted products chilled against the ground, inhibiting zeolite crystallization.
During rapid deposition of the thick, wet succession thermal insulation allowed the persistence of elevated temperatures for
a time sufficient for enhancement of hydration-dissolution processes in the volcanic glass. The highly reactive alkali-trachytic
glass quickly buffered the acid pH of the system, favouring phillipsite crystallization followed by chabazite nucleation.
The variable zeolite content reflects fluctuating emplacement conditions (e.g. changes in water content and temperature).
Cooling of the upper and relatively thin distal deposits inhibited the zeolitization process, thereby preserving the primary
unlithified deposit.
Received: 25 May 1999 / Accepted: 28 October 1999 相似文献
The American crocodile was declared endangered in the United States in 1975. At that time 75% of the remaining crocodile nests were in Everglades National Park, in Florida Bay. In 1980, the National Park Service established a crocodile sanctuary in northeastern Florida Bay to protect nesting and nursery habitat. In 1985, a monitoring program, focused on nesting, growth, and survival, was established to evaluate the effects of modified water deliveries on crocodiles in Florida Bay. The number and range of crocodile nests increased between 1970 and 1995, but nesting success decreased slightly. Nests on artificial substrates in the Greater Flamingo-Cape Sable area accounted for most of the increase in nests. Nests on artificial substrates were more prone to predation by raccoons. At least 1.5% of marked hatchlings survived for more than 12 mo, and growth rates were variable. Detailed information on growth and survival of crocodiles is still lacking. It is no longer a question of whether crocodiles with survive in Florida Bay, but how ecosystem restoration and management can be applied to improve conditions for crocodiles. 相似文献
Stochastic optimization methods, such as genetic algorithms, search for the global minimum of the misfit function within a given parameter range and do not require any calculation of the gradients of the misfit surfaces. More importantly, these methods collect a series of models and associated likelihoods that can be used to estimate the posterior probability distribution. However, because genetic algorithms are not a Markov chain Monte Carlo method, the direct use of the genetic‐algorithm‐sampled models and their associated likelihoods produce a biased estimation of the posterior probability distribution. In contrast, Markov chain Monte Carlo methods, such as the Metropolis–Hastings and Gibbs sampler, provide accurate posterior probability distributions but at considerable computational cost. In this paper, we use a hybrid method that combines the speed of a genetic algorithm to find an optimal solution and the accuracy of a Gibbs sampler to obtain a reliable estimation of the posterior probability distributions. First, we test this method on an analytical function and show that the genetic algorithm method cannot recover the true probability distributions and that it tends to underestimate the true uncertainties. Conversely, combining the genetic algorithm optimization with a Gibbs sampler step enables us to recover the true posterior probability distributions. Then, we demonstrate the applicability of this hybrid method by performing one‐dimensional elastic full‐waveform inversions on synthetic and field data. We also discuss how an appropriate genetic algorithm implementation is essential to attenuate the “genetic drift” effect and to maximize the exploration of the model space. In fact, a wide and efficient exploration of the model space is important not only to avoid entrapment in local minima during the genetic algorithm optimization but also to ensure a reliable estimation of the posterior probability distributions in the subsequent Gibbs sampler step. 相似文献
The Italian “Guidelines for the seismic risk classification of constructions” approved in February 2017 define the technical principles for exploiting tax deductions with respect to seismic strengthening interventions on existing buildings (Sismabonus). Tax deductions represent a unique opportunity to improve the seismic safety of the existing Italian building stock. The guidelines are very simple and allow practitioners to deal with the sophisticated concepts behind modern seismic design, such as expected annual losses (EAL) and repair costs (expressed as a fraction of the Reconstruction Cost: %RC). The seismic risk classes of buildings and the class upgrade due to strengthening interventions can be assessed using the principles included in the guidelines. The seismic risk class is the minimum between the class defined by the building safety index at the ultimate limit state and the one related to the EAL. The latter class depends on the area under the curve of the expected losses, which is easily obtained by computing the safety index converted in the return period (annual frequency) at different limit states and the relevant %RC. This paper illustrates the technical principles at the base of the guidelines and the procedure used to calibrate the repair costs associated with the different limit states using the actual repair costs monitored in the reconstruction process following recent Italian earthquakes. Finally, simple tools to estimate the cost of the strengthening interventions to improve the seismic capacity at the life-safety limit states are provided. 相似文献
The fronts of tunnels excavated under particularly difficult ground conditions by employing conventional tunnelling methods are commonly supported: the stabilization is usually achieved either by improving the mechanical properties of the soil (injections, jet grouting, soil freezing, etc.) or by introducing linear inclusions. This last technique, consisting in the introduction of pipes (usually made of fibreglass reinforced polymers) in the front, is particularly popular since it is very simple to adapt the reinforcement geometry, length and number to the different conditions encountered during the excavation. The design of this reinforcement technique is nowadays based on very simplified approaches: on either empirical formula or the limit equilibrium method. In a previous paper, the authors numerically studied the mechanical response of unreinforced fronts in cohesive soils and defined a non-dimensional front characteristic curve. In this paper, the authors intend to take into consideration the role of reinforcements by following the same approach. A procedure allowing the definition of the reinforced non-dimensional front characteristic curve, once the reinforcement pattern is assigned, is introduced. The practical use of this curve is described.
Batter piles are widely used in geotechnical engineering when substantial lateral resistance is needed or to avoid the interference with existing underground constructions. Nevertheless, there is a lack of fast numerical tools for nonlinear soil‐structure interactions problems for this type of foundation. A novel hypoplastic macroelement is proposed, able to reproduce the nonlinear response of a single batter pile in sand under monotonic and cyclic static loadings. The behavior of batter piles (15°, 30°, and 45°) is first numerically investigated using 3D finite element modeling and compared with the behavior of vertical piles. It is shown that their response mainly depends on the pile inclination and the loading direction. Then, starting from the macroelement for single vertical piles in sand by Li et al (Acta Geotechnica, 11(2):373‐390, 2016), an extension is proposed to take into account the pile inclination introducing simple analytical equations in the expression describing the failure surface. 3D finite element numerical models are adopted to validate the macroelement that is proven able to reproduce the nonlinear behavior in terms of global quantities (forces‐displacements) and to significantly reduce the necessary computational time. 相似文献
At sub‐arc depths, the release of carbon from subducting slab lithologies is mostly controlled by fluid released by devolatilization reactions such as dehydration of antigorite (Atg‐) serpentinite to prograde peridotite. Here we investigate carbonate–silicate rocks hosted in Atg‐serpentinite and prograde chlorite (Chl‐) harzburgite in the Milagrosa and Almirez ultramafic massifs of the palaeo‐subducted Nevado‐Filábride Complex (NFC, Betic Cordillera, S. Spain). These massifs provide a unique opportunity to study the stability of carbonate during subduction metamorphism at P–T conditions before and after the dehydration of Atg‐serpentinite in a warm subduction setting. In the Milagrosa massif, carbonate–silicate rocks occur as lenses of Ti‐clinohumite–diopside–calcite marbles, diopside–dolomite marbles and antigorite–diopside–dolomite rocks hosted in clinopyroxene‐bearing Atg‐serpentinite. In Almirez, carbonate–silicate rocks are hosted in Chl‐harzburgite and show a high‐grade assemblage composed of olivine, Ti‐clinohumite, diopside, chlorite, dolomite, calcite, Cr‐bearing magnetite, pentlandite and rare aragonite inclusions. These NFC carbonate–silicate rocks have variable CaO and CO2 contents at nearly constant Mg/Si ratio and high Ni and Cr contents, indicating that their protoliths were variable mixtures of serpentine and Ca‐carbonate (i.e., ophicarbonates). Thermodynamic modelling shows that the carbonate–silicate rocks attained peak metamorphic conditions similar to those of their host serpentinite (Milagrosa massif; 550–600°C and 1.0–1.4 GPa) and Chl‐harzburgite (Almirez massif; 1.7–1.9 GPa and 680°C). Microstructures, mineral chemistry and phase relations indicate that the hybrid carbonate–silicate bulk rock compositions formed before prograde metamorphism, likely during seawater hydrothermal alteration, and subsequently underwent subduction metamorphism. In the CaO–MgO–SiO2 ternary, these processes resulted in a compositional variability of NFC serpentinite‐hosted carbonate–silicate rocks along the serpentine‐calcite mixing trend, similar to that observed in serpentinite‐hosted carbonate‐rocks in other palaeo‐subducted metamorphic terranes. Thermodynamic modelling using classical models of binary H2O–CO2 fluids shows that the compositional variability along this binary determines the temperature of the main devolatilization reactions, the fluid composition and the mineral assemblages of reaction products during prograde subduction metamorphism. Thermodynamic modelling considering electrolytic fluids reveals that H2O and molecular CO2 are the main fluid species and charged carbon‐bearing species occur only in minor amounts in equilibrium with carbonate–silicate rocks in warm subduction settings. Consequently, accounting for electrolytic fluids at these conditions slightly increases the solubility of carbon in the fluids compared with predictions by classical binary H2O–CO2 fluids, but does not affect the topology of phase relations in serpentinite‐hosted carbonate‐rocks. Phase relations, mineral composition and assemblages of Milagrosa and Almirez (meta)‐serpentinite‐hosted carbonate–silicate rocks are consistent with local equilibrium between an infiltrating fluid and the bulk rock composition and indicate a limited role of infiltration‐driven decarbonation. Our study shows natural evidence for the preservation of carbonates in serpentinite‐hosted carbonate–silicate rocks beyond the Atg‐serpentinite breakdown at sub‐arc depths, demonstrating that carbon can be recycled into the deep mantle. 相似文献